Theory of coherent acoustic phonons in InGaN/GaN mult i- quantum wells 
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A microscopic theory for the generation and propagation of 
coherent LA phonons in pseudomorphically strained wurtzite 
(0001) InGaN/GaN multi-quantum well (MQW) pin diodes 
is presented. The generation of coherent LA phonons is driven 
by photoexcitation of electron-hole pairs by an ultrafast Gaus- 
sian pump laser and is treated theoretically using the density 
matrix formalism. We use realistic wurtzite bandstructures 
taking valence-band mixing and strain-induced piezoelectric 
fields into account. In addition, the many-body Coulomb in- 
teraction is treated in the screened time-dependent Hartree- 
Fock approximation. We find that under typical experimen- 
tal conditions, our microscopic theory can be simplified and 
mapped onto a loaded string problem which can be easily 
solved. 

PACS Number(s): 63.10.+a, 63.20.Kr, 63.22.+m 



I. INTRODUCTION 

In recent years, experiments have shown that optical 
excitation of electron-hole pairs in semiconductors by ul- 
trafast lasers can coherently excite, longitudinal optical 
phonon modes in scmiconductorsJj113 In uniform bulk 
semiconductors, since the laser wavelength is much larger 
than the lattice spacing, the photo-generated carriers are 
typically excited by the optical pump over spatial areas 
that are much larger than the lattice unit cell. As a 
result, the excited carrier populations are generated in 
a macroscopic state and the carrier density matrix has 
only a q ss Fourier component. Coupling of the pho- 
toexcited carriers to the phonons leads only to coherent 
optical phonon modes with q w 0. Since the frequency 
of the {SiO acoustic phonon is zero, coherent acoustic 
phonons are not excited in bulk semiconductors. 

In semiconductor superlattices, even though the laser 
pump has a wavelength large compared to the lattice 
spacing, the pump can preferentially generate electron- 
hole pairs in the wells. The result is to create photoex- 
cited carrier distributions that have the periodicity of 
the superlattice. Since the density matrix of the photo- 
excited carrier populations now has ag^O Fourier com- 
ponent, the photo-excited carriers can not only couple 
to the optical phonon modes, but they can also gen- 
erate coherent acoustic phonon modes with a nonzero 



frequency and wavevector q w 2n/L where L is the su- 
perlattice period. In superlattices, the coherent phonon 
oscillation of zone folded acoustic phonons has been ob- 
served in AlAs/GaAs superlatticesBt£l. However, the 
reflection modulation, observed to he on the order of 
AR/R ~ 1CT 5 - 1CT 6 , is very.small.Ea 

Recently, C.-K. Sun et. alia reported studies of co- 
herent acoustic phonon oscillations in wurtzite (0001) 
InGaN/GaN multi-quantum well samples with strain 
induced piezoelectric fields. Owing to the strong piezo- 
electric fields at the interfaces, huge coherent acoustic 
phonon oscillations were observed. The oscillations were 
strong enough to be seen in the transmision (rather than 
the usual reflectivity) with AT/T ~ 10~ 2 - 10~ 3 . The 
oscillation frequency, in the THz range, corresponding to 
the LA phonon frequency with q ss 2tt/L, varied between 
samples in accordance with their different superlattice 
periods, L. 

In this paper, we formulate a microscopic model for 
the generation of coherent acoustic phonons in strained 
wurtzite superlattices via ultrafast laser photo-excitation 
of real carriers. Whereas in bulk systems the micro- 
scopic theory of coherent LO,-phonons can be mapped 
onto a forced oscillator modelcl, we show that coherent 
LA phonon generation in superlattices, under appropri- 
ate conditions, can be mapped onto a loaded string model 
which is readily solved for the lattice displacement. Since 
acoustic phonons are almost the same in the well as in the 
barrier, _to lowest order we can treat the string as being 
uniformllj. The forcing term on the string, however, is 
not uniform since photoexcitation of carriers occurs only 
in the wells. 

Our paper thus provides justification for using a sim- 
ple, uniform string model with a nonuniform forcing 
term, rather than a more complicated microscopic theory. 
In addition, we provide a microscopic expression for the 
forcing term to use in the simplified string model. The 
string model provides additional insight into the physics 
of the coherent LA phonons. 



II. MICROSCOPIC THEORY 

In this section, we derive the microscopic theory for 
coherent acoustic phonon generation in superlattices 
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and multiple quantum wells, including the effects of (i) 
bandstructure, (ii) strain, (iii) piezoelectric fields, (iv) 
Coulomb interactions, and (v) laser optical excitation. 
In section III, we will show how this reduces to a sim- 
plified driven uniform string model with a nonuniform 
forcing term and a microscopic expression for the forcing 
term. 

We model photogeneration of electrons and holes and 
the subsequant excitation of coherent acoustic phonons in 
a multi-quantum well (MQW) pin diode shown schemati- 
cally in Fig. [I]. The intrinsic active region consists of a left 
GaN buffer region, several pseudomorphically strained 
(0001) In x Ga\^ x N quantum wells sandwiched between 
GaN barriers, and a right GaN buffer region as indi- 
cated in the figure. The P and N regions are assumed to 
be abruptly terminated p- and n-doped GaN bulk layers 
separated by a distance, L, across which a voltage drop, 
Ay = Va, is maintained. Photoexcitation of carriers 
is achieved by means of an ultrafast laser pulse incident 
normally along the (0001) growth direction, taken to co- 
incide with the z-axis. 



\v, 2) = ^L\{X-tY)])-^\{X + iY) |) (3b) 
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\vA) = ~\(X + iY) T) iY) I) (3d) 
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The phase factors, a and /?, are functions of the angle 
4> = tan -1 (k y / k x ) and are given by 



a(0) 



V2 



o i(37r/4+30/2) 



(4a) 



A. Bulk Bandstructure 

In bulk systems, the conduction and valence bands 
in wurzite crystals including the effects of strain are 
treated using effective mass theory. Near the band edge, 
the effective mass Hamiltonian for electrons is described 
by a 2 x 2 matrix which depends explicitly on electron 
wavevector, k, and the strain tensor, e. The electron 
Bloch basis states are taken to be 



\c,l) = \S]) 



|c,2) = |Si>. 



(la) 



(lb) 



The conduction band Hamiltonian is diagonal and we 
have (relative to the bottom of the conduction band) 
(Refs. H and fil) 



2;„2 



m 



(M = { 



n 2 k? 

T 9m* 
;)} 12x2- 



2m* z 



(2) 



where 12x2 is the identity matrix. The electron effective 
masses along z (taken parallel to the c-axis) and in the 
xy plane are m* and m* yl respectively, k\ = k x + ky 7 and 
£xx, tyy and e zz are strain tensor components, and a C)Z 
and a CtXy are the deformation potentials. 

The Hamiltonian for the valence bands is a 6x 6 matrix. 
Following Ref. [l5|, the Hamiltonian (relative to the top 
of the valence band) can be block diagonalized into two 
degnerate 3x3 submatrices if we adopt the Bloch basis 
states 

\v, 1) = -^=\(X +iY)\) + ^=\(X- iY) |) (3a) 
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J(tt/4+4,/2) 



(4b) 



The block diagonalized Hamiltonian can be written as 
H^siKe) 



(5) 



ff 3 L x 3 (k,e 
where the upper and lower blocks of the Hamiltonian are 



and 



^(M) 



F K t -iH t 
K t G A - iH t 
iH t A + iH t A 



F K t iH t 
K t G A + iH t 
-iH t A -iH t A 



(6a) 



(6b) 



The elements appearing in the 3x3 Hamiltonian ma- 
trices are 



F = Ai + A 2 + A 



G = Ai - A 2 + A 



2m, 



-A 5 k 2 t 



H = A 6 k t k z 

ZrriQ 



A = V2A 3 



(7a) 
(7b) 

(7c) 

(7d) 
(7e) 
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{A x kl + A 2 k\) + D 1 e zz + D 2 ( ) (7f) 



(A 3 k 2 z + A 4 k 2 ) + D 3 e zz + D 4 {e xx + e yy ) (7g) 



InEq. (|),the^s are effective mass parameters, the D[s 
are the Bir-Pikus deformation potentials, and the A's are 
related to the crystal field splitting, A cr , and spin-orbit 
splitting, A so , by Ai = A cr and A 2 = A 3 = A so /3. m 
is the free electron mass. 



B. Quantized carrier states in MQW diodes 

In quantum confined systems such as the pin diode 
shown in Fig. |l|, we must modify the bulk Hamiltonian. 
The finite MQW structure breaks translational symme- 
try along the z direction but not in the xy plane. Thus, 
quantum confinement of carriers in the MQW active re- 
gion gives rise to a set of two-dimensional subbands. The 
wavefunctions in the envelope function approximation 
are 



C k (r) 



E 
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\oi,J) 



(8) 



where a = {c, v} refers to conduction or valence sub- 
bands, n is the subband index, k = {k x , fc y , 0) = (fc, 4>) is 
the two-dimensional wavevector, and j labels the spinor 
component. For conduction subbands, (a=c) j = 1,2 
while for valence subbands (a=v) j = 1...6. The slowly 
varying envelope functions F" k ^(z) are real and depend 
only on fc = |k|, while the rapidly- varying Bloch basis 
states | a, j) defined in Eqs. (|l|) and (j3|) depend on <fi in the 
case of valence subbands as given in Eq. (Q). The area of 
the MQW sample in the xy plane is A, and p = (x, y, 0) 
is the projection of r in the plane. 

The envelope functions satisfy a set of effective-mass 
Schrodinger equations 

{Hfj, (fc) + hr Wa (*) - K (k)] } K tkd , (z) = o, 



jj' 



subject to the boundary conditions 



o) = K,kA z 



L) = 0. 



(9) 



(10) 



where L is again the length of the MQW diode structure 
(c.f. Fig. [l]), V a (z) are the quantum confinement poten- 
tials for conduction and valence electrons. E"(k) are the 
energy eigenvalues for the nth conduction or valence sub- 
band. Note that in the envelope function approximation, 
the subband energy depends only on the magnitude fc of 
the transverse wavevector and not on the angle <j). For 
the quantum confined case, the matrix operators H" -, (fc) 



depend explicitly on z and are obtained by making the re- 



placement k z 



■ dz 



and letting all material parameters 



be z-dependent operators in the matrices H a (k, e) given 
in Eqs. (^|) and (^|). To ensure the Hermitian property-of 
the Hamiltonian, we make the operator replacementst3 



d 







d 



dz' 



and 



B[z) 



dz 



B(z) 



8_ d_ 

dz dz 



B(z) 



(11a) 



(lib) 



The quantum confinement potentials V c {z) and V v (z) 
arise from (i) bandgap discontinuities between well and 
barrier regions, (ii) the strain-induced piezoelectric field, 
and (iii) the time-dependent electric field due to photoex- 
cited electrons and holes. Thus, 

V a (z,t) = V^gap^) + Vpiezo(z) + Vphoto(M), (12) 

and the bandstructure is explicitly time-dependent. 

Material parameters for InN and GaN used in this 
work can be found in Tabic |. For In x Ga\^ x N alloys, 
we interpolate between the GaN and InN values listed 
in the table. 

We obtain electron effective masses by linearly inter- 
polating the reciprocals of the masses as a function of the 
indium concentration x, i.e., the concentration dependent 
effective masses are taken to be 



1 



1 



m %y\ x ) \ m *xyJ InN 

and 



+ (1-*) 



1 



/)r; 



X V' GaN 



1 



m* z [x) 



z / InN 



(13a) 



(13b) 



z / GaN 



For the alloy band gap, E g (x), we use an expression 
incorporating a bowing paremeter: 

E g {x) = xE gJnN + (1 - x)E g ^ GaN - b x(l - x) , (l4) 

where the bowing parameter, b = 1.0 eT^.@ 

For all other material parameters, we use linear inter- 
polation in x to obtain values for the alloy. Since we can't 
find deformation potentials for InN, we use GaN values 
by default. In the absence of values of for either GaN 
or InN, we use linear interpolation in x to obtain eo and 
simply take eoo ~ eo- 

If the z-dependent band gap in the MQW is E g (z) as 
determined from Eq. ( |l4| ) and the Indium concentration 
profile, and E g ^ m { n — m\n z [E g (z)\ is the minimum band 
gap in the structure, then the confinement potentials for 
conduction and valence electrons are defined as 

Vc, E&P {z) = B 9 ,min + Qc {E g (z) - Eg tmin ) (15a) 



K,gap(z) = -(1 - Qc) (E g {z) - -E^min) 



(15b) 
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where the conduction band offset is taken as Q c = 0.6.E2I 
With these definitions for the V a , gap (z), the zero of the 
gap confinement potential is placed at the top of the va- 
lence band profile. 

The confinement potentials due to the strain-induced 
piezoelectric field are given by 



V piezo (z) = -\e\ E° z (z), 



(16) 



where |e| is the electric charge and E z (z) is the strain- 
induced piezoelectric field. In a pseudomorphically 
strained MQW diode, the bulk source and drain (as- 
sumed to have identical composition) are unstrained 
while the in-plane MQW lattice constants adjust to the 
source and drain values. For a MQW grown along [0001] 
(the z-direction) the z-dependent strain isE3 



,{z) = e vy {z) 



<20 — 0,(z) 

a(z) 



(17) 



Here ao is the lattice constant in the source and drain 
and a(z) is the z-dependent lattice constant iirfie MQW. 
Minimizing the overall strain energy, we findE3 



C 33 (z) ' xx{Zh 



(18) 



where C\ 3 (z) and 633(2:) are z-dependent elastic con- 
stants. 

There are several issues concerning strain that one can 
worry about. One is the critical well thickness beyond 
which the strain relaxes. Studies have showr£3 that with 
10 % Indium, and 6 nm well width, a 350 kV/cm field 
is measured, implying that the well is not fully relaxed 
and justifies using the pseudomorphic strain approxima- 
tion as we do in this paper for wells having 6 % Indium 
and thickness near 4 nm. For thick wells, this pseudor- 
morphic strain approximation clearly will begin to break 
down and a more detailed model will be needed. In- 
terface roughness can also play a role. The roughness 
will not significantly affect the acoustic phonon modes 
(see introduction) but may affect the photogeneration of 
carriersE-3. For simplicity, we do not consider interface 
roughness. 

The strain-induced polarization directed along the z- 
direction is given by 

P°(z) = e 31 {z) (e xx (z) + e yy {z)) + e 33 (z) e zz (z), (19) 

where e 3 i(z) and 633(2) are z-dependent piezoelectric 
constants. The unscreened piezoelectric field in the diode 
is obtained from the requirement that the electric dis- 
placement vanishescj. Thus, 



47T 



(P°(z) + P ), 



(20) 



where Pq is a constant polarization induced by exter- 
nally applied voltages and £0(2) is the position-dependent 
static dielectric constant. The value of Pq is obtained 



from the voltage drop across the diode (of length L) in 
the unscreened limit, i.e. with no photoexcited carriers. 
In this limit, the voltage drop between source and drain 
due to the induced piezoelectric field is just 



V A = ~ I dz E° z (z), 
10 



(21) 



from which Pq can be determined. 

When photoexcited electrons and holes are generated 
by the laser, then there is an additional time-dependent 
confinement potential, 



V photo (z,t) = -\e\ Ef oto (z,t). 



(22) 



This potential is obtained by solving the Poisson equa- 
tion in the diode for £?P hoto (z, t) subject to the boundary 
condition 



V A = - I dz E z otal {z,t) 



(23) 



Here, E z ° tal (z,t) is the total electric field and is just the 
sum of the strain induced electric field and the field due 
to photogenerated electrons and holes, i.e., 



E 



total 



(z,t)^E° z (z) + E^°(z,t). 



(24) 



Finally, we can write an effective-mass Schrodinger 
equation for the conduction electron envelope functions 
in terms of an effedctive electron potential, V^*' (z): 



8 



1 



8 



2 [dz m*(z) dz J '"'^ 
+ {K e// M ~ K(k)} K,kA z ) = °» (25) 
where the effective electron potential is 

VS"{z) =V c (z) + -^ f - 
2 m* y (z) 

+ a c .z(z)e zz (z) + a c ^ xy (z) (e 

xx (-2*) ~t~ £yy 

(z)) . (26) 

Similar expressions can be derived for valence electrons. 
We arrive at a set of coupled ordinary differential equa- 
tions (ODE's) subject to the two-point boundary value 
condition of Eq. (|l0|) . These are solved for the envelope 
functions and subband energies. 

In practice, we introduce a uniform grid, {zi}, along 
the z-direction and finite-difference the effective mass 
Schrodinger equations to obtain a matrix eigenvalue 
problem which can be solved using standard matrix 
eigenvalue routines. The resulting eigenvalues are the 
subband energies, E%(k), and the corresponding eigen- 
vectors are the envelope functions, F" k j(zi), defined on 
the finite difference mesh. 
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C. Second quantized electron hamiltonians 

We next describe the second quantized Hamiltonians 
for electrons moving freely in the MQW interacting via 
a screened Coulomb potential. We denote creation and 
destruction operators for electrons in conduction and va- 
lence subbands by k and c ai7l) k, respectively. The 
second quantized Hamiltonian for free electrons and holes 
is simply 



7~teQ — 



a,n,k 



k*--a,n,k- 



(27) 



The Coulomb interaction Hamiltonian is given by 

*- = 5£ E E 

a,n,k a' ,n' k^O 
X C Ln,k-K C l' : n',k'+K Ca ' ,n' M> Ca,n,k- (28) 

Eq. ( pq ) describes two-body interactions where electrons 
in states |a, ra, k) and |a', n', k') scatter to subband states 
|a,n,k — k) and |a',ra',k' + k), respectively. Note that 
to simplify things, we have neglected terms correspond- 
ing to Coulomb induced interband transitions (the "di- 
agonal approximation" ) since these are very unfavorable 
energetically.E£l The electrons thus stay in their original 
subbands (though they may scatter off of other electrons 
in different subbands) and exchange crystal momentum 
k. The matrix elements describing the strength of these 
transitions are given by 



Va.n, k / \ 
a',n',k'( K ) 



dz J dz 1 Vj re | [z — z 1 ) 

^2 F n,\k- K \,j( Z ) F n,k,j( Z ) 



(29) 



From Eq. (|2^), it is apparent that the symmetry relation 

(30) 



V a >n>M>( K ) = V a.«, k ( ) 



must hold. 

The Fourier transform in the xy plane of the screened 
Coulomb potential depends only on k = |k| and \z\ and 
is given by 



K(z) 



27re 2 er K ^ 



e A k e s (n) 



(31) 



To describe screening, we adopt an effective pseudody- 
namic dielectric function of the form 



1 



c s (k) k + k s 



(32) 



In the pseudodynamic screening model, we completely 
neglect screening by the massive holes and treat screening 



by the lighter conduction electrons in the static screening 
limit. The screening wavevector, k s , is computed in the 
2D limit. ThusB 



27re 2 dN' 



2D 



CO 



dfj, 



(33) 



where N2D, the two-dimensional conduction electron 
density, is related to an effective chemical potential, /z, 
by 



m* k B T ^ 



1 + exp 



k R T 



(34) 



In Eq. (p4[), E^(0) is the conduction subband energy eval- 
uated at k = 0. In our simulation, the value of fi is ob- 
tained by requiring that N2D, evaluated using Eq. (|34|), 
be equal to the value 



AMf)=E/«( M ) 



(35) 



obtained from the time dependent conduction electron 
distribution functions, /^(fc,i). 



D. Photogeneration of carriers 

Electron-hole pairs are created by the pump laser and 
we treat the electric field of the laser in the semiclas- 
sical dipole approximation. In this approximation, the 
electron-laser interaction Hamiltonian is 



H eL = -\e\ E 



(t). £ [d- (k)<i 



\^v^n' ,k H~ h.C. 



n.n' . k 



(36) 



where h.c. denotes the Hermitian conjugate of the first 
term. The laser field is E(t) and the dipole matrix ele- 
ments are 



<?„<(k) E D ;::- ;oi / dz w*) p n> 

3,3' 



uA z ) ■ (37) 



The vector operator, D^'J, (</>), is a 2 x 6 matrix with x, 
y, and z components. Thus, 

D$, (0) = D°* (</>) x + 1^(0) y + D^(</>) z , (38) 

where x, y, and z are unit vectors and 

P2 



i P2 



a -(3* a -(3* 
-a* (3 a* -(3 



-a -(3* -a -(3* 
-a* -(3 a* (3 



(39a) 



, (39b) 
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DT(<f>) 



E n 



0-/300 (3 
-0* -P* 



(39c) 



The 6x2 vector operator, Dp c . (</>), is related to the 2x6 
operator, D^'",(^), by 



(40) 



In Eq. (|3^) , a and /3 are the ^-dependent phase factors 
defined in Eq. (^) and the Kane parameters, P\ and P 2 , 
for wurtzite materials are related to the effective masses 
and energy gaps byc3 



m 



- 1 



2mo \m 
{E g + A 1 +A 2 ){E g + 2A 2 )-2Aj 
E q + 2A 2 



(41a) 



P? = 



mo 



- 1 



2m \m* y 
E g {(E g + Ax + A 2 ){E g + 2A 2 ) - 2A|} 
(S g + Ai + A 2 )(S g + A 2 )-A2 



(41b) 



For the semiclassical laser field, we write the real elec- 
tric field as 



£(t) 



(42) 



where to is the photon frequency, e is a complex unit 
polarization vector, and £{t) is the pulse shape envelope 
function. We assume a Gaussian pulse shape 



£(t) = £q exp 



t-t 



1 

2 In 2 . 



(43) 



centered at t = to with an intensity full width at half 
maximum (FWHM) of r. The maximum electric field 
strength, £$, is related to the pump fluence, J- ', by 



£0 = 



16tt T /In 2 



c n^T 



(44) 



where n w is the index of refraction at the photon fre- 
quency. 

For linearly polarized light incident normally on the 
MQW, the polarization vectors are real and given by ei- 
ther x or y. For circularly polarized light, the polariza- 
tion vectors are complex and given byESI 



x ± iy 
V2 



(45) 



In Eq. (flq), the upper sign refers to left circularly polar- 
ized light (positive helicity) and the lower sign refers to 
right circularly polarized light (negative helicity). 



E. Coupling to LA phonons 

We treat the acoustic phonons in the MQW as bulk- 
like plane-wave states with wavevector q. Since the sys- 
tem exhibits cylindrical symmetry, only q = q z lon- 
gitudinal acoustic phonons are coupled by the electron- 
phonon interaction. The free LA phonon Hamiltonian 
can be written as 



Ti-Ao = hu q b\b q 



(46) 



where 6j and b q are creation and destruction operators for 
LA phonons with wavevector q = q z. The wave vector 
component, q, of LA phonons in the MQW is thus defined 
in an extended zone scheme where — 00 < q < 00. The 
phonon dispersion relation is given by a linear relation 



C s \q\ 



' C33 



(47) 



where po is the mass density and C s is just the LA phonon 
sound speed for propagation parallel to zB3 In comput- 
ing the LA sound speed in the linear phonon dispersion 
relation of Eq. (E^), we neglect the z-dependence of the 
material parmeters and use bulk GaN values for C33 and 
Pa- 

The LA phonons in wurtzite MQW's interact with 
the electrons through deformation potential and screened 
piezoelectric scattering. The electron-LA phonon inter- 
action in an MQW is governed by the Hamiltonian 

(48) 

This Hamiltonian describes the scattering of an electron 
from subband state |a,n',k) to subband state \ot, n, k) 
with either the emission or absorption of an LA phonon. 
We note that the electron wave vector, k, in the xy plane 
is conserved in this process since, as noted earlier, the 
phonon wave vector in the xy plane is zero. 

The interaction matrix elements describing deforma- 
tion and screened piezoelectric scattering are 



i( l v n,n'{k,q) 



2p (twj q ) V 

47r|e| e 33 



(49) 



where V is the crystal volume. The first term in Eq. 
( p9| ) desribes deformation potential scattering while the 
second term describes screened piezoelectric scattering. 

The relative strengths of the various transitions are de- 
termined by form factors for deformation potential and 
piezoelectric scattering. The form factor for screened 
piezoelectric scattering is given by 
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K,n>(k, q) = E f dZ KkA Z ) ^ F n'^A Z ) > ( 5 °) 
3 J 

while the form factor for deformation potential scattering 
is defined to be 

j J 

(51) 

The form factor for deformation potential scattering 
is similar to the form factor for piezoelectric scattering 
except that in summing over spinor components, j, the 
terms are weighted by j-dependent deformation poten- 
tials 6" which can be representd by the row vectors 

@j = { a c , z , a Cj2 } , (52a) 

9} = {D l + D 3 ,D 1 + D 3 , D u D l + D 3 , D l + D 3 , D x } . 

(52b) 

for conduction and valence electrons, respectively. 

F. Electron density matrices 

We define statistical operators in terms of the electron 
and phonon eigenstates. The electron density matrix is 

<(M) = (cl,„ !k (t)c a :^ k (t)), (53) 

where ( ) denotes the statistial average of the non- 
cqulilbrium state of the system. 

The interband components of the density matrix, 
^nn' 4 ) t) an d -^n'nO 4 ' describe the coherence be- 
tween conduction and valence electrons in subbands n 
and n' and are related to the optical polarization. The 
intraband components of the density matrix, N (k, i) 
describe correlations between different subbands of the 
same carrier type if n ^ n'. If n = n', N^'"(k,t) = 
/™(k, t) is just the carrier distribution function for elec- 
trons in the subband state, "0"k( r )' defined in Eq. (|J). 

G. Coherent phonon amplitude 

The coherent phonon amplitude of the q-th phonon 
mode, \q), is defined to beQ 

D q (t) = ( b q (t)+bl q (t) ) . (54) 

The coherent phonon amplitude is related to the macro- 
scopic lattice displacement, U(z,t), and velocity, V(z,t), 
through the relations 



U ^=^2^V eiq ' D ^ (55) 

V(z, t) = V J f , ^ (56) 

The coherent phonon amplitude, D q (t), will vanish if 
there are a definite number of phonons in the mode, i.e., 
if the phonon oscillator is in one of its energy eigenstates, 
\q). In this case, there is no macroscopic displacement of 
the lattice. _ 

The coherent phonon distribution isQ 

N<°\t) = (6j(t)) ( bq (t)) (57) 

and the total phonon distribution, J\f q (t), can be sepa- 
rated into coherent and incoherent contributions as fol- 
lows; 

M q {t)= (b q (t) bl q {t) ) 

= N™ h {t)+N™ coh {t). (58) 

In general, a mode can have a number of both coherent 
and incoherent phonons, but only the coherent phonons 
contribute to the macroscopic lattice displacement. 

We note that at the beginning of the experiment, there 
are no coherent phonons present, i.e. Af q oh (t) = 0, and 
the incoherent phonon population is described by a ther- 
mal distribution, A/"™ co,l (i) ~ e -^ 9 /k B T - 

H. Equations of motion 

In this section we develop equations of motion for the 
electron density matrices and coherent phonon ampli- 
tudes. The electron density matrices obey the general 
equations of motion 

^ = < \H, < n , k c QVl McJ ) , (59) 

where [ ] denotes the commutator and ( ) denotes the av- 
erage over an initial ensemble. The density matrices are 
defined in the electron picture and initially the valence 
bands are filled while the conduction bands are empty. 
We have f°(k,t = -co) = and f%(k,t = -oo) = 1 
which implies 

N n,n> ( k > 1 = _ °°) = S n,n> S a>v $a> ,v ■ (60) 

The total Hamiltonian, 7i, is the sum of the Hamiltonians 
described in the previous sections, i.e. 

7~L = tleO + W ee + 'HeL + ~HaO + 7~teA ■ (61) 

In deriving equations of motion for the density matri- 
ces, we make the ansatz that the density matrices depend 
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only on k = |k|. We use the rotating wave approximation 
(RWA) to factor out the rapid e lut behavior of the inter- 
band density matrix elements, N^ v n ,{k,t). In the RWA, 
we have 



(62) 



where N% n >(k,t) is a slowly varying envelope function. 
In addition, we treat the Coulomb interaction in the 
time-dependent Hartree-Fock approximation by factor- 
ing four-operator averages arising from TL e L into appro- 
priate products of two-operator averages as described in 
Ref . |7|. 

The resulting equations of motion for the density ma- 
trices are 

m 

+ JE' {K m (k)N c m %,(k) NZ%(k)A c m>n ,(k)} (63a) 

m 

dNn:n dt {k - = l - s v Ak)} K : :,(k) 

m 

+ { A l m (k)N v m : n ,(k) - NA m {k)Ka..n'(k)} (63b) 



A" / (fc,t) = E« ,(k,t)+Q° ,(k,t) 



(65) 



The first term in the renormalization energy matrix ( |65| ) 
is the generalized exchange self-energy matrix arising 
from the Coulomb interaction and is given by 



,(M) 



E 

k'^k 



V 



,(|k-k'| 



(66) 



where V"'™/ k ,( |k — k'| ) are angular averaged Coulomb 
interaction matrix elements. The second term in Eq. ( |65| ) 
accounts for renormalization due to coupling of carriers 
to coherent acoustic phonons. We have 



?n,n'(k> t) 



Y,D q (t) M° n ,(k,q), 



(67) 



where D q (t) is the coherent phonon amplitude and the 
electron-phonon matrix elements, M" n ,(k,q), are de- 
fined in Eq. (|4|) . The self energy corrections in Eq. ( |6^ ) 
are small, though the can be important in some circum- 
stances. 

In computing the angular averaged Coulomb matrix el- 
ements in Eq. (p9), we assume small momentum transfer, 
k, and use the fact that the envelope functions, k j(z), 
depend weakly on k to obtain an effective interaction, 



Va.nl / \ 
a>',fc'( K ) 



E 



F n,(k,k'),j 



dz I dz' V K (z — z') 



AW 



I7>a 

r n',(k,k') 



(68) 



9N ^ k,t) = l - {£ c n {k) - £l,{k) - hu}N^,{k) 

-iJ2 {a^(k)N^ n ,(k) - NAAmZAk)} 

rn 

+ JE' '{ A n,Ak)N C A n ,(k)~ NAAkWra.Ak)} • (63c) 

m 

The equations of motion for N^ ,c ^,(k,t) are redundant 

since ^ / (fc,t) = (JV^yfc,t))*." 

The first terms on the right-hand side of Eq. (|6|) de- 
scribe the free oscillation of the density matrices in the 
renormalized singel particle energy bands. The time- 
dependent single particle energies are 



*E(M)= W) + A" (M) 



(64) 



where E"(k) are the single particle subband energies 
in the absence of conduction electrons and holes and 
A" n (k,t) describes the time-dependent renormalization 
of the single particle subbands. 

The renormalization energies, A% n (k,t), are the di- 
agonal elements of a generalized renormalization energy 
matrix (in the subband indices) 



where, by definition, 

t n,(k,k'),j\ Z )\ 



(69) 



The effective Coulomb interaction, V, defined in Eq. (|6£ 
is an even function of k and is symmetric in a and k, thus 
preserving the symmetry relation (30). Preserving this 
symmetry is essential in order to maintain conservation 
of carriers in the scattering process. 

The second terms in Eq. ( |63| ) describe photoexcitation 
of electron-hole pairs by the pump laser. The system 
reacts to an effective field which is the sum of the applied 
field and the dipole field of the electron- hole excitations. 
This gives rise to a matrix of generalized Rabi frequencies 
in the subband indices 



E 



muk) 



V£&(|k- 



£{t) 



d c AAk) + 



(70) 



k'\)NAAk',t), 



which can be shown to satisfy the symmetry relations 

™z>(k) = (nnyjk))* (7i) 
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The Gaussian pump envelope function, £(t), is defined 
in Eq. (k3|) and the optical dipole matrix elements 



(k) (72) 



are angular averages in the xy plane of the vector dipole 
matrices dotted into the polarization vector. From 
Eq. ([37]), the ^-dependence of d^ v n , (k) only appears in 
Dj'J/ (</>) and we can get the angular averages by setting 
a(4>) = a avg = (1 - i)/3it and /3(0) = /3 Q1 , S = (1 + i)/ 71 " m 
the 2 x 6 matrices D^ v '(</>), D%. v (<j>) and D c /(4>) defined 
in Eq. @. 

The last terms in Eq. ( |63| ) are similar in structure to 
the renormalization corrections in the Hartree-Fock ener- 
gies but are more complicated due to mixing among sub- 
bands and involve the off-diagonal components of A" n , . 
The prime on the summation sign indicates that terms 
containing factors of N"'", (k) are excluded from the sum 
since these terms have already been incorporated in the 
renormalized Hartree-Fock energies in Eq. (|64|). 

The coherent phonon amplitudes, D q (t), satisfy the 
driven harmonic oscillator equations 



d 2 D q {t) 
dt 2 



)DS) = 



2uj n 



E 



{<£(M)-<W, S nX }, (73) 



subject to the initial conditions 



Dg(t = "CO) = 



dD q (t = -CO) 

dt 



= 0. 



(74) 



The closed set of coupled partial differential equations, 
( |63| ) and ( 73 ) , for the carrier density matrices and coher- 
ent phonon amplitudes are converted into a set of coupled 
ordinary differential equations (ODE's) by discretizing k 

and q and solving for AT"'", (ki) and D{qi) for each of 
the mesh points fcj and qi . The resulting initial value 
ODE problem is then solved wsing a standard adaptive 
step-size Runge Kutta routinetj. 

The phonon distributions do not appear in the coupled 
set of equations (^) and (|73]), If necessary, they can be 
determined from N (fc) and the pair of equations 



g^fcoh 



dt 



|lm J2 MlAk,q)B q N^,{k) (75a) 

a,n,n' ,k 



and 



dB q 
dt 



q^q 



(75b) 



In Eq. ( 75t ) , B q (i) = (b q (t)} satisfies the initial condition 
B q (t = —00) — 0. For the incoherent phonon distribu- 
tion, 



QJ^fincoh 



(76) 



so no incoherent phonons are generated and the inco- 
herent phonon population maintains its initial thermal 
equilibrium distribution. 



III. LOADED STRING MODEL 

The microscopic equations are rather daunting and de- 
tailed. In this section, we show how they can be simpli- 
fied (under certain conditions) to a more tractable model, 
manely that of a driven uniform string, provided one uses 
the appropriate driving function, S(z, t), which is nonuni- 
form. The microscopies, including details of the super- 
lattice band structure and photogeneration process are 
included within the driving function. 

In our detailed numerical simulations, we use the full 
microscopic formalism discussed in the previous sections. 
However, we gain a lot of insight if we can deal with the 
lattice displacement, U(z,t), directly. If we assume that 
the acoustic phonon dispersion relation is linear as in Eq. 
(p7|) , then we find that U(z, t) satisfies the loaded string 
equation 



d 2 U(z,t) 



ct 



d 2 U(z,t) 



dt 2 s dz 2 

subject to the initial conditions 



S(z,t) 



dU (z, t = —00) 

U(z, t = -00) = ^— = 0. 

v ' dt 



(77) 



(78) 



The LA sound speed, C s , is defined in Eq. (|47|), and the 
driving function, S(z,t), is given by 



I2h C s \q\ 



a,n,n k,<? i 

x [N%Z(k,t) - 8 a , v 6 n>n ,} 



(79) 



One may question whether a linear phonon dispersion 
relation is valid in a superlattice. For small wavevector, 
q, for which elasticity theory holds, the dispersion rela- 
tion for LA phonons in a superlattice is linear with a 
dispersion, u> = C s q, where C s is the "averaseT sound 
speed of LA phonons in the well and barrierstHl. This, 
in fact, has been experimentally verified in InGaNVGaN 
superlattice samples studied by C.-K. Sun et. alO. 

Note that coherent acoustic phonon generation in a 
superlattice is qualitatively different than coherent opti- 
cal phonon generation in a bulk system where only the 
gsO optic mode can be excited. As a result, both the 
amplitude, U(z,t), and the Fourier transform of the am- 
plitude, D q (t), for an optic mode in bulk satisfy a forced 
oscillator equation. For the nonuniform, multiple quan- 
tum well case, one can excite acoustic modes with q =/= 0. 
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The Fourier transform of the amplitude D q (t) of a coher- 
ent acoustic phonon obeys a forced oscillator equation, 
but owing to the linear dependence of u)(q) on q ,the am- 
plitude itself, U(z,t), obeys a 1-D wave equation with a 
forcing term, S(z,t). 

Another important point is that Eq. ( fn\) can be taken 
to be a uniform string with a non-uniform forcing func- 
tion. This is because the speed of sound is approximately 
the same in both the GaN and InGaN layers (a more 
detailed theory would take into account differences in 
the sound velocities in each layer). For propagation of 
acoustic modes one can neglect, to lowest order, the dif- 
ferences between the different layers (this is not true for 
the optic modes). The non- uniformity of the forcing func- 
tion S(z, t) results from differences in the absorption (not 
sound velocity) in the well and barrier layers and is there- 
fore z dependent. We thus see from Eq. (77) that un- 
derstanding coherent acoustic phonons in multiple quan- 
tum wells is equivalent to understanding a uniform string 
with an inhomogeneous forcing term, S(z,t), containing 
the microscopies. 

To simplify Eq. ([79]), we neglect valence band mixing 
and assume the effective masses, sound speeds, and cou- 
pling constants are uniform over regions where S(z, t) ^ 
0, i.e., in regions where carriers are being photogener- 
ated. We also assume the pump pulses are weak enough 
so that screening of the piezoelectric interaction can be 
neglected. Finally, if the pump duration is long enough 
so that transient effects associated with photogeneration 
of virtual carriers can be ignored, then the off-diagonal 
elements of the carrier density matrices in Eq. ([79]) can 
be dropped. In this case, the driving function takes on 
the simple form 



SCM) = X>CM), 



(80) 



where the summation index, v, runs over carrier species, 
i.e., conduction electrons, heavy holes, light holes, and 
crystal field split holes. 

Eq. ( |8"(i| ) says that each carrier species makes a separate 
contribution to the driving function. The partial driving 
functions, S v (z,t), are 

S v (z, t) = ±- [a® + 47r|e| 633 } p u (z, f), (81) 

PO I OZ Eoo J 

where the plus sign is used for conduction electrons and 
the minus sign is used for holes. Here p v (z, t) is the pho- 
tognerated electron or hole number density, which is real 
and positive, and po is the mass density. We note that 
the loaded string equation for the propagation of coherent 
phonons together with the simplified driving function in 
Eqs. ( p0| ) and ( |8l| ) have also been independently derived 
by other authors in the limit 633 = O.EH 

In Eq. (|8"l|), the partial driving function for a given 
species is obtained by applying a simple operator to the 
photogenerated carrier density. This operator is a sum 



of two terms, the first due to deformation potential scat- 
tering and the second to piezoelectric scattering. The 
piezoelectric coupling constant, 633, is the same for all 
carrier species, while the deformation potential, a„, de- 
pends on the species. For conduction electrons, a v = a c , z , 
for heavy or light holes, a v = D\ + D3, and for crystal 
field split holes, a v = D%. 

It is interesting to note that Planck's constant does not 
appear in either the loaded string equation, ( [77]), or in 
its associated driving function defined in Eqs. (pQ) and 
(pl|). Thus, we find that coherent LA phonon oscillations 
in MQW's can be viewed as an essentially classical phe- 
nomenon, an observation that was made in the context of 
coherent LO phonon oscillations in bulk semiconductors 
by Kuznetsov and Stanton in Ref. |] 

The driving function, S(z,t), satisfies the sum rule 



dz S(z,t) = 0. 



(82) 



This is most easily seen from Eqs. (|80|) and (|8l|). but it 
also holds for the general expression in Eq. Q79|). The 
significance of the sum rule is readily appreciated. After 
the pump dies away, the carrier density in Eq. (81), ne- 
glecting tunneling between wells, is essentially constant 
and thus S(z, t) is time-independent. In the loaded string 
analogy, the integral of the driving function over position 
is proportional to the average force per unit length on the 
string. If this integral were non-zero, then the center of 
mass of the string would undergo a constant acceleration 
resulting in the buildup of an infinite amount of kinetic 
energy. Such an alarming result in the context of co- 
herent LA phonons is precluded by the sum rule in Eq. 

For a given driving function, the wave equation (|77|), 
together with the initial conditions ([78]), can be solved 
for the coherent phonon lattice displacement using the 
Green's function method.tEl Thus, 



OO pQO 



OO J — OC 



U{z,t)= dt' dz' G(z- z',t-t') S(z',t'). 



(83) 

In our MQW diode model, the substrate is assumed to 
be infinite and the Green's function in this case is just 



G(z,t) 



2 a 



{&(z + C s t)-e(z-C a t)}, (84) 



where Q(x) is the Heaviside step function. 

We note that the loaded string model described above 
is not restricted to the special case of an infinite sub- 
strate and can be extended to study the generation and 
propagation of coherent LA phonons in more complicated 
heterostructures. If the driving function, S(z,t), due to 
photoexcited carriers is localized, then the assumptions 
leading to Eqs. ( |80| ) and ( pi] ) need only hold in those re- 
gions where S(z, t) in nonvanishing. The wave equation 
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applies in regions where the LA sound speed, C s , is con- 
stant. Heterostructure in which the LA sound speed is 
piecewise constant have abrupt acoustic impedance mis- 
matches which can be handled by introducing more com- 
plicated Gpea's functions or by using other standard 
techniques .Etlj An example of such a problem would be a 
MQW structure embedded in a free standing substrate in 
which coherent LA phonons generated in the MQW could 
bounce back and forth between two parallel substrate-air 
interfaces. 



IV. RESULTS 

In this section, we discuss simulations based on our 
microscopic theory of coherent LA phonon generation in 
a pin diode structure with four periods of InGaN/GaN 
MQW's photoexcited by a Gaussian pump normally in- 
cident along the (0001) z-direction. The parameters for 
our numerical example are listed in Table g The MQW 
dimensions and Gaussian pump parameters were chosen 
to match those typically encountered in room tempera- 
ture pump-probe differential transmission measurements 
of coherent iA phonon oscillations carried out by C.-K. 
Sun, et. alliil on In om Ga .viN / GaN MQW structures 
having fourteen periods. 

A. Bulk wurtzite bandstructure 

Bulk wurtzite GaN and InN are direct gap materi- 
als with band gaps of 3.4 and 1.95 eV, respectively. The 
bulk bandstructure of unstrained wurtzite GaN is shown 
in Fig. ^. As can be seen from equations (||) and (^), 
the bandstructure is anisotropic and depends on k z , the 
wavevector along the (0001) z-axis, and k t , the wavevec- 
tor within the xy plane perpendicular to the z-axis. The 
effective mass conduction band is two-fold degenerate 
and has a parabolic dispersion with anisotropic effective 
masses m* = 0.19 along the z-direction and m* y = 0.18 
in the xy plane. 

The two-fold degenerate valence bands are mixtures of 
heavy hole (HH) , light hole (LH) , and crystal-field splitoff 
hole (CH) character. At the zone center, the off-diagonal 
components of the 3x3 upper and lower Hamiltonians 
in Eq. (^|) vanish and the valence bands can be labeled 
according to their pure state wavefunctions at k = 0. For 
the zone-center HH state, the degenerate wavefunctions 
are the basis states \v, 1) and \v,4) defined in Eq. (||). 
For the zone-center LH state, the wavefunctions are \v, 2) 
and \v,b), and for the CH band the zone-center wave- 
functions are \v,3) and \v,6). The heavy-hole effective 
masses along z and xy are m^ H = \Ai + A^ 1 = 1.96 
and m^ y H = \A 2 + A 4 - A 5 | _1 = 1.92 for heavy holes, 
m% H = mf ff = 1.96 and = \A 2 +A i +A^\- 1 = 0.14 
for light holes, and m^ H = = 0.14 and = 

1^42 1 — 1 = 1-96 for crystal-field splitoff holes. 



B. Pseudomorphic strain 

Bulk GaN and InN have different lattice constants so 
when an (0001) InGaN MQW structure is grown, a sig- 
nificant lattice mismatch occurs between the In x Ga\^ x N 
wells and GaN barriers. For the InGaN MQW diode 
specified in Table [n[ we assume pseudomorphic strain 
conditions. In a pseudomorphically strained device, the 
lattice constant throughout the MQW adjusts to the 
value of the lattice constant in the bulk N and P sub- 
strates in order to minimize the overall strain energy. In 
our simulated diode, the substrates are n- and p- doped 
GaN, so the lattice constant throughout the device takes 
on the GaN value, i.e. ao = 3.189 A. The non-vanishing 
position-dependent strain tensor components, e xx , e vv , 
and e zz , for the MQW diode, as computed from Eqs. ( |l7| ) 
and (|l8|), are shown in Fig. |as a function of z. Clearly, 
the GaN barriers are unstrained since the N and P sub- 
strates are composed of GaN and all the strain from the 
lattice mismatch is accomodated in the In0.0QGa0.Q4N 
wells. 



C. Built-in piezoelectric field 

The presence of strain in the MQW's results in the 
creation of a strain-induced polarization, P®{z), directed 
along z as described by Eq. ([l9|). The strain-induced 
polarization, in turn, results in a strong bult-in piezo- 
electric field which can be computed from Eqs. ( p0| ) and 
(^l|), given the strain field and the dc bias, Va, applied 
across the diode. The computed strain-induced piezo- 
electric field, E®(z), and the piezoelectric confinement 
potential, V p i czo (z), which result from the strain field in 
Fig. H are shown in Fig. [|. Prior to the application of the 
pump pulse, we assume the applied dc bias, Va, has been 
adjusted so flat-band biasing in the diode is achieved, i.e. 
Va is such that the band edges seen in Fig. |^ are periodic 
functions of position. 

Given the piezoelectric field and confinement poten- 
tials, position dependent band edges for the MQW can 
be computed. The conduction and valence band edges for 
our pseudomorphically strained MQW diode are shown 
as functions of position in Fig. [|. These are just the con- 
finement potentials, V a (z) — V r Qiga p(z) + Vpi czo (z), in the 
diode prior to photoexcitation. It is clear from Fig. |E| 
that the confinement of elecrons and holes in the MQW 
is mostly due to strong built-in piezoelectric fields which 
result in the triangular confinement potentials seen in 
each well. 



D. Photogeneration of carriers 

In our numerical example, we simulate photoexcitation 
of electrons and holes and the generation and subsequent 
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propagation of coherent LA phonons in the hypotheti- 
cal MQW diode when a Gaussian pump laser pulse is 
normally incicent along the z-axis. As seen in Table || 
the Gaussian pump pulse is assumed to be left circularly 
polarized with a photon energy of 3.21 eV. The pump 
fluence is taken to be 100.0 pJ/cm 2 and the Gaussian 
FWHM is taken to be 180.0 fs. The expermient is as- 
sumed to take place at room temperature. 

In Fig. ^, the computed conduction and valence sub- 
band energies are shown as functions of k for the 
In x Ga\- x N diode. At the chosen pump energy of 
3.21 eV, electrons from the first two valence subbands 
are excited into the lowest lying conduction subband. 

The computed densities of photoexcited electrons and 
holes, neglecting and including Coulomb interaction ef- 
fects, are shown as functions of position and time in Figs, 
fj] and ||, respectively and the total photoexcited electron 
density per unit area as a function of time is shown in Fig. 
||. In Fig. ^| the pulse shape is shown for comparison. We 
find that including Coulomb effects decreases the total 
photogenerated carrier density. The electrons and holes 
screen the built-in piezoelectric field widening the effec- 
tive band gap. This quantum confined Stark effect acts 
to suppress the photogeneration of carriers. 



E. Generation of coherent phonons 

The driving function, S(z,t), for the driven string 
equation (0) is shown in Fig. |l^ as a function of po- 
sition and time. The driving function has units of accel- 
eration and in Fig. [l^, we compute S(z,t) using the full 
microscopic formalism of Eq. ( 79 ) . 

For comparison, we also computed the driving func- 
tion in the simplified loaded string model of Eqs. (pQ) 
and ( |8l| ) using the carrier densities shown in Fig. |i|to 
facilitate the comparison. Since the photoexcited holes 
are predominantly a mixture of heavy- and light- holes, 
we use a v = D\ + D 3 in computing hole deformation po- 
tential contributions in Eq. (pl|). The sum over species, 
v, then yields the total driving function 



b{Z,t)= — < a c , z — H ) pelec(Z,t) 

PO I OZ £oo J 

--((£>i + D 3 ) I- + Pbole ( z>t ), (85) 

P0 I OZ Coo J 

where p c i cc (z, t) and phoie(z, t) are the total conduction 
electron and valence hole densities plotted in Fig. |l The 
resulting S(z,t) is shown in Fig. |ll|. 

By comparing Figs. [l0|and|ll|, we see that for the diode 
structure and Gaussian pump used in our simulation the 
simplified loaded string model produces essentially the 
same results as those obtained using the full microscopic 
formalism. 

Acoustic LA phonon generation due to the piezoelec- 
tric effect depends on the piezoelectric constant, 633, the 



number of photogenerated electrons and holes, as well 
as the spatial separation of electron and hole densities 
brought about by the strong built-in piezoelectric field in 
the MQW's. From Eq. (p5[), the piezoelectric contribu- 
tion to the driving function is given by 

Spiczo^, t) = — {p c lcc(z, t) - Phdle(z, t)} (86) 

P0 e oo 

In the absence of a built-in piezoelectric field (such as 
found in a square well with infinite barriers) we would 
have /9 e lec(M) ~ Phoie(z,t) and hence S piczo (z,t) w 
even for relatively large values of 633. The built-in 
piezoelectric field serves to spatially separate the elec- 
trons and holes so that p B \ ec (z,t) 7^ Phoio(z,t) and hence 
<S'piezo(^, t) 7^ 0. However, if the built-in piezoelectric field 
is too strong and the spatial separation of electrons and 
holes too large, then p e \e C (z,t) — Pho\ c {z 7 t) « 0. This is 
because the overlap between the conduction and valence 
envelope functions enters into the optical dipole matrix 
elements in Eq. j3^) . If there's negligable overlap be- 
tween electron and hole envelope functions due to strong 
piezoelectric fields then d™ n , (k) sa 0, no electron-hole 
pairs are photogenerated and once again S p i czo (z^ t) « 0. 

The deformation potential contribution to the driving 
function is given by 



dp c \ cc {z, t) (Di + D 3 ) dp ho i c (z, t) 



Po 



dz 



Po 



dz 



(87) 



From Table ||, the conduction electron deformation po- 
tential, a c 2, is roughly twice the valence hole deforma- 
tion potential, Di + D 3 . Thus, the two terms in Eq. fls^) 
are of comparable magnitude. The first term, due to con- 
duction electrons, gives rise to a contribution to Sd e {(z, t) 
which is localized on the right side of each MQW while 
the second, due to valence holes, gives rise to contribution 
which is localized on the left hand side of each MQW. 

In our simulation, we find that piezoelectric and defor- 
mation potential contriburtions to the driving function 
are comparable. This is seen in Fig. [l^ where S P i CZO (z) 
and Sdef(z), along with their sum, are plotted at t = 2 ps. 
In this example, we find that Sp lezo (z) makes the domi- 



nant contribution to S(z,t) as can be seen in Fig. 12 



The macroscopic lattice displacement, U(z, t), and ve- 
locity field, V(z,t), can be obtained from the coherent 

Al- 



phonon amplitudes, D q (t), using Eqs. (55) and ( p6[ ) 
ternatively, they can also be obtained from the driving 
function, S(z,t), through the Green's function solution 
of the driven string equation. In Fig. O, we plot the 
macroscopic lattice displacement, U(z, t), and velocity 
field, V(z,t), for coherent LA phonon modes generated 
by the driving function shown in Fig. fTol 
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F. Coherent phonon energy 

From the lattice displacement, U(z,t), we can obtain 
the total energy density per unit volume associated with 
coherent LA phonons, 

S LA (z, t) = T LA (z, t) + V LA {z, t), (88a) 

as the sum of a kinetic energy density term, 

p (z) fdu(z,ty 2 



T ^M- 2 V dt 
and a potential energy density term, 



V LA (z,t) 



C 33 (z) fO U(z,t) 
dz 



(88b) 



3c) 



The coherent LA phonon energy per unit area, E^ A (t), 
is obtained by integrating £L A (z,t) over position, z: 



E LA {t) 



dz £ LA (z,t). 



(89) 



The total energy density, £L A (z,t), for coherent LA 
phonons as a function of position and time is shown in 
Fig. |IJ and in Fig. |l5| the total energy density is plotted 
as a function of position for equally spaced values of the 
time ranging from t = to t = 8 ps in increments of 2 ps. 
The curves in Fig. [l5| for different times have been offset 
to avoid overlapping. At short times, the evolution of the 
total phonon energy density is complicated, but the long 
times behavior, t > 6 ps, can be easily understood. 

As t — > oo, a localized energy density appears in the 
MQW region due alm ost entirely to the potential en- 
ergy term in Eq. (88c). This is due to near steady-state 
loading by the driving function at long times. Assuming 
the driving function, S(z,t), is approximately constant 
at long times, the loaded string equation, (^7[), can be 
integrated once in the steady state limit. We find the 
steady state solution 

dU(z) ^ 
dz ~ 



dz 



CI ' 



(90) 



from which the long time behavior of the LA phonon 
energy density per unit volume in the MQW's, 



C 33 (z) 
2 A 



dz 



c? 



(91) 



can be obtained. The fact that the energy density in Eq. 
( |9l| ) is localized in the MQW's follows directly from the 
sum rule ( |S2| ) and is clearly seen in Figs. [L4| and [| 

In addition to the localized energy density, which re- 
mains behind in the MQW's, two propagating wave 
trains consisting of four pulses each are seen to exit 
the MQW region and travel off to infinity at the acous- 
tic phonon sound speed, C s . The distance between the 



pulses is just the inter-well separation distance. In these 
radiating wave trains, the kinetic and potential energy 
densities, Tl A (z, t) and VL A {z,t), are found to be equal 
as one would expect. 

The power spectrum of the coherent LA phonon energy 
density in g-space can be written in terms of the coherent 
phonon amplitudes D q (t). The power spectrum for the 
total coherent LA phonon energy density, 



£ LA {q,t) = T LA (q,t) + V LA (q,t), 
is again the sum of a kinetic energy term, 



dD q (t) 



Of 



and a potential energy term 
1 



V LA {q,t) 



2A 



hWq\D q (t)\' 



(92a) 



(92b) 



(92c) 



The phonon energy density per unit area is obtained 
by summing the power spectrum over positive phonon 
wavevectors, q. Thus, 



E LA (t)=J2 £L A (q,t). 

q>0 



(93) 



The total energy density power spectrum for coherent LA 
phonons as a function of phonon wavevector, q, and time 
is shown in Fig. |l6|. The peak near q = is associated 
with buildup of the steady state energy density localized 
in the MQW region. Secondary peaks are seen near qo = 
0.59 nm~ 1 and twice this wavevector, i.e. q\ — 2 qo — 
1.18 nm . The wavevector qo corresponds to the the 
wavevector of the MQW period, 



2tt 



qo = 



Li, 



(94) 



where L w and are the well and barrier widths. 

The total coherent LA phonon energy pe r unit area can 
be obtained from either Eqs. ( |89| ) or (p3|). In Fig. [Tt], 
we show the total phonon energy per unit area, E^ A (t), 
as a function of time for the coherent LA phonons gener- 
ated by the driving function shown in Fig. [l(]. The total 
energy per unit area is the sum of kinetic- and potential- 
energy terms. For comparison, the pulse shape is shown 
as a dotted line. It is clear from the figure that the 
buildup of energy in coherent LA phonons takes place on 
a time scale that is much longer than the pump duration. 
In addition, we see that the total energy buildup in the 
phonons saturates at around 5 ps and that some strong 
but rapidly decreasing oscillations are superimposed on 
top of an increasing trend. 

The saturation phenomenon results from the fact that 
we have a finite number of quantum wells and not an 
infinite superlattice. The results can best be explained 
in terms of the driven string equation. In general, the 
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rate at which energy is fed into the phonon. field per unit 
area is described by the energy equation,E3 



dE hA 
dt 



Po 



dz S(z,t) V(z,t), 



(95) 



in which S(z, t) and V(z, t) are the phonon driving func- 
tion and velocity fields defined in Eqs. (|79| ) and ([56|), and 
Po is the GaN mass density used for computing the sound 
speed, C s , in Eq. j47|). The energy equation simply says 
that the rate at which energy is added to a driven string is 
proportional to the local force times velocity integrated 
over the length of the string. The integral in Eq. j95| ) 
vanishes when the transient velocity field, V(z,t), exits 
the MQW region in which the driving function, S(z,t), 
is localized. Thus, the time, i S at, required for El A to 
saturate is just the time it takes for an LA sound wave 
to cross the MQW, i.e. t s&t ~ W/C s , where W is the 
width of the MQW region over which the driving func- 
tion is localized. In our example the LA sound speed is 
C s = 80 A/ps in GaN and the MQW width (four well 
and three barrier layers) is W = 38lA, from which we 
obtain i sat = 4.8 ps. 

The oscillations of Ela observed in Fig. |l^ reflect the 
number and periodicity of the diode MQW's. The pump 
laser generates spatially periodic electron and hole distri- 
butions, as seen in Fig. |7J, due to the fact that the pump 
photoexcites carriers in the wells but not the barriers. 
From each of the wells in the MQW, two sound pulses 
emerge traveling in opposite directions thus giving rise to 
an outwardly propagating velocity field pattern, V(z,t), 
with four peaks traveling outward in each direction as 
seen in Fig. 13. The driving function, S(z,t), on the 
other hand, is localized in the MQW's and is relatively 
constant in time after the pump pulse dies out. The 
driving functions localized in each well do work on four 
traveling velocity disturbances, the one generated in the 
well itself as well as the ones generated in the three neigh- 
boring wells that subsequently pass by. This gives rise 
to the three peaks plus saturation plateau seen in Fig. 
p"7j . The time interval between coresponding peaks in ad- 
jacent wells is just the time it takes LA sound waves to 
travel between wells, i.e. £ pcr iod = {Liu + Lf,)/C s , where 
L w and L are the well and barrier thicknesses. For the 
MQW structure, L w = 63.0 A and L b = 43.0 A and 
we have i per iod = 1-325 ps which agrees with the peak- 
to-peak time between the first and second peaks seen in 
Fig. |i~7l 



V. SUMMARY AND CONCLUSIONS 



ence of strain in the MQW's results in the creation of 
strain-induced built-in piezoelectric fields on the order 
of sevaral MV/cm which significantly alter the electronic 
and optical properties of the diode structure. In partic- 
ular, the effective band gap can be lower than the band 
gap in unstrained InGaN wells due to the presence of 
triangular piezoelectric potentials. 

To a first approximation, the generation of coher- 
ent LA phonons is driven by optical photoexcitation of 
electron-hole pairs by an ultrafast Gaussian pump laser. 
Under typical experimental conditions, the propagation 
of coherent LA phonons is described by a uniform loaded 
string equation for the lattice displacement where the 
time- and position-dependent driving force on the string 
is a function of the photoexcited carrier density. This dif- 
fers from coherent LO phonon oscillations in bulk systems 
where the coherent LO phonons obey a forced oscillator 
equation. Both deformation potential and piezoelectric 
coupling mechanisms contribute to the driving force. We 
find that deformation potential coupling contributes a 
driving force proportional to the derivative of the carrier 
density while piezoelectric coupling contributes a driving 
force proportional to the photoexcited carrier density. 

We found that the driving term in the loaded string 
equation is suddenly turned on by rapid generation of 
electron-hole pairs by the pump and remains approxi- 
mately constant theafter. This sudden displaceive load- 
ing of the string results in a new static equilibrium lattice 
displacement. This new static equilibrium displacement 
corresponds to a population of coherent LA phonons with 
5»0. As the lattice adjusts to the new equilibrium, co- 
herent LA phonons are transmitted in the positive and 
negative z-directions at the LA sound speed. These trav- 
eling coherent LA phonons are characterized by q « 2tt/L 
where L is the superlattice period. 

The formalism described here can be applied to the 
analysis of more complicated device geometries as well 
as more complicated laser pulse sequences. This gives a 
simpler method for calculating the coherent LA phonon 
generation in more complicated geometries and gives ad- 
ditional insight into the acoustic coherent response. 
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We have developed a microscopic theory for the gener- 
ation and propagation of coherent LA phonons in pseudo- 
morphically strained wurzite (0001) InGaN/ GaN multi- 
quantum well (MQW) pin diodes. Both GaN and InN 
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TABLE I. Material parameters for wurtzite InN and GaN. 
Material parameters for In x Gai- x N are obtained through 
interpolation in x as described in the text. 



Parameter 



InN 



GaN 



Lattice constants 
a ( A) 
co( A) 

MO 

Direct band gaps (eV) 

E~a 



Electron effective masses (mo) 



Hole effective mass parameters 

Tx 

A 2 

A, 

A 4 

A s 

A 6 

Hole splitting energies (meV) 
Ai = A cr 
A 2 = A s0 /3 
A 3 

Electron deformation potentials (eV) 

0>c,xy 

a>c,z 

Hole deformation potentials (eV) 

Lh 

D 2 

D 3 

Piezoelectric constants (C/m 2 ) 

e3i 

e33 



Elastic stiffness constants (GPa) 
Cu 

Cl2 
C13 



3.540 
3.708 
0.377 



1.95 



0.10 c 
0.11 c 



-9.28 c 
-0.60 c 
8.68 c 
-4.34 c 
-4.32 c 
-6.08 c 



17.0 c 
1.0 c 
1.0 c 



-0.57 
0.97 



190 
104 
121 



3.189 
5.185 
0.376 



3.40 



0.18 c 
0.19 c 



-7.24 c 
-0.51 c 
6.73 c 
-3.36 c 
-3.35 c 
-4.72 c 



22.0 c 
3.67 c 
3.67 c 



0.7 ' 
2.1 ' 
1.4 ' 
-0.7 



-0.49 
0.73 



374 
106 
70 
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C33 
C44 

Static dielectric constant 

£0 



182 
10 ' 



15.3 



a Ref. 
b Ref. 
c Ref. 
d Ref. 
e Ref. 
f Ref. 
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TABLE II. Simulation parameters for photogeneration of 
coherent acoustic phonons in a four well MQW diode under 
flat band biasing conditions. A schematic of the diode struc- 
ture is shown in Fig. |l| . 

MQW diode structure 



Left GaN buffer width (A) 

Number of wells 
Well width (A) 
Indium fraction in well 
GaN barrier width (A) 

Right GaN buffer width (A) 

Applied bias 
Va {V) 

Lattice temperature 
YjK) 

Pump parameters 
Photon energy (eV) 
Fluence ( fiJ/cm 2 ) 
Gaussian FWHM (fs) 
Polarization 



43.0 

4 
63.0 
0.06 
43.0 

43.0 



-0.261 



300.0 



3.21 
160.0 
180.0 
Left circular 



FIG. 1. Schematic diagram of the In x Gai- 
turn well diode structure. 
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FIG. 2. Bulk GaN valence band structure using effective 
mass parameters taken from Table [l| The bands are plot- 
ted along the (0001) fc z -axis and along the transverse fc t -axis 
within the xy plane. The anisotropic zone-center effective 
masses for heavy holes (HH), light holes (LH), and crystal 
field splitoff holes (CH) are indicated. 
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FIG. 3. Strain tensor components for pseudomorphically 
strained In x Ga\- x N mult i- quantum well diode as a function 
of position. The diode parameters are listed in Table [□]. 
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FIG. 4. Electric field and potential for the strain field in 
Fig. [| The applied dc bias, Va, has been adjusted so 
flat-band biasing is achieved, i.e. so that the band edges 
are periodic functions of position. The diode parameters are 
listed in Table [l[ 
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FIG. 5. Conduction and valence band edges for pseudo- 
morphically strained In x Ga\~ x N multi-quantum well diode 
as a function of position. The applied dc bias, Va, has been 
adjusted so flat-band biasing is achieved, i.e. so that the band 
edges are periodic functions of position. The diode parame- 
ters are listed in Table [o]. 
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FIG. 6. Conduction and valence subband energies as func- 
tions of k for the In x Ga\- x N diode structure described in 
Table 
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(a) Electrons 




(b) Hole: 




FIG. 7. Density of excited carriers computed in the ab- 
sence of Coulomb effects for (a) electrons and (b) holes as 
functions of position for the In x Gai~ x N diode structure and 
laser pumping parameters shown in Table |l| 
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FIG. 8. Density of excited carriers including Coulomb ef- 
fects for (a) electrons and (b) holes as functions of position 
for the In x Ga\- x N diode structure and laser pumping pa- 
rameters taken from Table [n|. 
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FIG. 9. Total photoexcited electron density with and with- 
out Coulomb effects as a function of time for the In x Ga 1 - x N 
diode structure and laser pumping parameters listed in Table 
[n]. The pulse shape (arb. units) is shown for comparison. 



— 



18 



4 



2 
O 



— 6 



" v\/ MM 




CO 



2.0 



CL 



ID 
I 

o 



N 

00 



30 
20 

10 


■10 



= 2 ps 



Total 




Piezo 



Deform = 



■40 -20 20 40 
Position (nm) 



FIG. 10. Driving function, S(z,t), for the coherent LA 
phonon wave equation as a function of position and time for 
the In x Gai- x N diode structure and laser pumping parame- 
ters in Table is computed using the full microscopic 
expression of Eq. (179) . 



FIG. 12. Driving function, S(z,t), in the simplified loaded 
string model at t = 2 ps for the coherent LA phonon wave 
equation as a function of position for the In x Gai- x N diode 
structure and laser pumping parameters in Table |TJ. The 
total driving function, S(z,t), is the sum of piezoelectric and 
deformation potential contributions, S , p i C zo(z, i) and Sdd(z, t). 




FIG. 11. Driving function, S(z,t), in the simplified loaded 
string model for the coherent LA phonon wave equation as 
a function of position and time for the In x Gai- x N diode 
structure and laser pumping parameters in Table O. 
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FIG. 13. Lattice displacement, U(z,t), and velocity field, 
V(z,t), for coherent LA phonons generated by the driving 
function shown in Fig. 
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FIG. 14. Total energy density, £l,a{z, t), for coherent LA 
phonons as a function of position and time for the driving 
function shown in Fig. [jl]. The total integrated energy den- 
sity as a function of time is obtained by integrating over po- 
sition, z. 
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FIG. 15. Total energy density, £la{z,£), for coherent LA 
phonons as a function of position for several values of the 
time, t, for the driving function shown in Fig. [j^. The curves 
for different times have been offset to avoid overlapping. 
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FIG. 16. Total energy density power spectrum, £LA{q,t), 
for coherent LA phonons as a function of phonon wavevector, 
q, and time for the driving function Fig. [jJJ. The total in- 
tegrated energy density as a function of time is obtained by 
integrating over the phonon wavevector, q. 
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FIG. 17. Integrated energy density, ELA{t), as a function 
of time for coherent LA phonons generated by the driving 
function shown in Fig. hd. The total integrated energy den- 
sity is the sum of kinetic- and potential- energy terms. The 
pulse shape (arb. units) is shown for comparison. 
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